typeset using JPSJ.sty <ver.l.Ob> 



Frustrated Spin System 
in 6»-(BEDT-TTF) 2 RbZn(SCN) 4 



00 
On 
On 



in 



c3 



i 

c 

o 
o 



> 

00 

in 
o 

oo 

On 
03 



O 

o 



X 



Hitoshi Seou and Hidetoshi Fukuyama 

Department of Physics, Faculty of Science, University of Tokyo, Tokyo 113 

(Received ) 

The origin of the spin gap behavior in the low-temperature dimerized phase of #-(BEDT- 
TTF)2RbZn(SCN)4 has been theoretically studied based on the Hartree-Fock approximation for 
the on-site Coulomb interaction at absolute zero. Calculations show that, in the parameter region 
considered to be relevant to this compound, antiferromagnetic ordering is stabilized between 
dimers consisting of pairs of molecules coupled with the largest transfer integral. Based on this 
result an effective localized spin 1/2 model is constructed which indicates the existence of the 
frustration among spins. This frustration may result in the formation of spin gap. 

KEYWORDS: BEDT-TTF, spin gap, Hartree-Fock approximation, on-site Coulomb interaction, antiferromagnetic 
spin ordering, frustration 



Two-dimensional (2D) organic conductors composed 
of BEDT-TTF (abbreviated as ET) have been attract- 
ing spread attention due to the variety of their electronic 
properties. In spite of the fact that the electron concen- 
tration of the conducting ET layer is fixed at 3/4 in com- 
pounds expressed as (ET^X, where X are anions, vari- 
ous phases such as paramagnetic metal, superconductor, 
paramagnetic insulator, antiferromagnetic (AF) insula- 
tor and even a spin-gapped (SG) insulator emerge. To 
understand these diverse behaviors, the role of Coulomb 
interaction is crucial because its value is considered to 
be of the order of the bandwidths in these compounds. 
Such effects of on-site Coulomb interaction have been in- 
vestigated systematically by Kino and Fukuyama within 
the Hartree-Fock (HF) approximation for k-(ET)2-X", a- 
(ET) 2 I 3 and a-(ET) 2 MHg(SCN) 4 BB ) It was shown by 
them that the ground state properties of ^-(ETphX can 
be explained by considering each dimer as a unitEP as was 
argued by Kanoda on experimental grounds.© The term 
'dimer' denotes the pair of parallel molecules between 
whom the transfer integral is noticeably larger than the 
other transfer integrals. Hence the AF insulating phase 
of k-(ET)2^ is considered to be a Mott insulator due to 
the existence of one hole per dimer. a If 

Such theoretical procedures of combing the HF calcu- 
lations and the quantum spin model so derived also allow 
us to understand the possible origin of the SG behavior 
in A-(BETS)2GaX^Y4_ z , whose donor is a seleno-analog 
molecule of ET, i.e. BETS (abbreviation for BEDT- 
TSeF) .ET The existence of a SG in this system had orig- 
inally been_suggested from magnetic susceptibility mea- 
surementSjEF and its origin was not due to the spin-Peierls 
(SP) transition. This compound has a dimeric structure 
similar to r-(ET)2A; the intradimer transfer integral is 
about twice that of the other transfer integrals. Thus 
the insulating phase can be considered as a 2D localized 
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spin- 1/2 system where each hole is localized on a dimer 
but an intrinsic alternation of the superexchange cou- 
plings between tbsse spins exists, which may lead to the 
formation of SGS 1 

Recently another 2D organic compound 0-(ET)2.RbZn(SCN)4 
has been synthesized which shows a SG behaviorJZP The 
existence oFa, SG has beeiLSuggested from magnetic sus- 
ccptibilityoLp and NMREf measurements. The 0-type 
ET compounds in general have a donor structure with 
two ET molecules in a unit cell as shown in Fig. 1(a) 
and have a simple 3/4-filled band, that is, they do not 
have a dimeric structure. However, in the case of 6- 
(ET)2RbZn(SCN)4, when it is slowly cooled, a metal- 
insulator transition .occurs at 190 K accompanied by 
a lattice modulation© which leads to the dimerization 
along the c-direction, thus the unit cell volume doubles 
and the transfer integrals become modulated as in Fig. 
1 (b) . No SP lattice distortion has been observed so farEr 
in the temperature range where the magnetic suscepti- 
bility shows a steep decrease. It can be seen in Fig. 
1(b) that molecule pairs (1-3) and (2-4) with intradimer 
transfer integrals t p z may be considered as dimers, but 
the value of t p s is not so large as compared to the other 
transfer integrals. Since in organic conductors such as 
(TMTSF) 2 X and k-(ET) 2 A the dimers consist of paral- 
lel molecules, we may even be able to consider the par- 
allel (1-2) and (3-4) molecule pairs as dimers possessing 
an intradimer transfer integral t c ±. 

The aim of this study is to elucidate whether dimers 
can be considered as units in the donor layer, and if it 
is so, which pairs of molecules should be considered as 
dimers. Furthermore, based on these results, the origin 
of the SG phenomenon will be investigated. 

In order to theoretically study the electronic struc- 
ture of the slowly cooled low temperature phase of 8- 
(ET) 2 RbZn(SCN) 4 , only its 2D donor plane (Fig. 1(b)) 
is considered where the transfer integrals between ET 
molecules are assumed as in Fig. 1(b) together with the 
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on-site Coulomb interaction U on ET molecule. The 
Hamiltonian is given as follows; 

H = ^ (tijah !' + h ' c ) + U ^ n i] n ili (!) 

<i, i 7> cr i 

where < i,j > denotes the neighboring site pair, a is 
the spin index which takes \ and J,, m a and a ia (aw) 
denote the number operator and the creation (annihila- 
tion) operator for the electron of spin a at the ith site, 
respectively. 

The Coulomb interaction U is treated in the HF ap- 
proximation in a manner similar to that in refs. 1 and 2. 
Our calculations are carried out at T = and the average 
electron density is fixed at 3/4- filling, namely one hole 
per two ET molecules as in the real compound. Several 
types of AF solutions are taken into account as shown 
in Fig. 2 and their energies are compared as a function 
of U so that the true ground state is searched for. The 
total energy £ is calculated as 

^ cell ,, 

Lka l 

where N ce u is the total number of cells, eika is the Ztli 
eigenvalue at each k and uf (e) is the Fermi distribution 
function. Due to the structure of intcrdimcr transfer 
integrals which we will discuss later, 3 different types 
of AF (a),(b) and (c), as shown in Fig. 2 exist, where 
the spins within the pairs (1-3) and (2-4) are parallel 
and AF ordering occurs between the pairs. In the case 
of (d) type AF ordering, the spins within the (1-2) and 
(3-4) pairs are parallel and 2D AF ordering takes place 
between pairs. 

The calculated energies of each state, a, relative to the 
paramagnetic solution, AE(a) = E a — E p , for four dif- 
ferent types of AF solutions are shown in Fig. 4. These 
results show that (a) and (b) types of AF, which are 
degenerate, are stabilized within realistic values of U, 
namely U = 0.6 — 1.0 eV. In these (a) and (b) types AF, 
the hole density on each ET molecule is similar but dif- 
fers by a few percent from +0.5, and the magnitude of 
spin moments is 0.8 — 0.9/is per dimer. Hence we can 
consider that the electrons are almost fully localized on 
each dimer. 

From the above results, we shall consider the pairs 
(1-3) and (2-4) as dimers. As shown in refs. 2 and 5 
the resulting dimer model can be constructed by esti- 
mating the superexchange couplings using the relation 
J ~ At 2 /V dimer {t is the effective interdimer transfer 
energy and V dimer is the effective Coulomb interaction 
between two holes in a dimer), which is given in Fig. 3 
with the magnitudes of J2 ~ 0.21Ji and J3 ~ 0.19Ji. 
This estimate shows that J\ is much larger than J2 and 
J3 is in accordance with the HF results where (a) or (b) 
type AF is stabilized. 

The quantum fluctuations among spins are taken into 
account in this model. The phase diagram of the ground 
state of such a 2D spin-1/2 model on the,plane of J2/J1 
and J3/J1, can be deduced as in Fig. 5.EJ When J 3 = 
(J2 = 0), the model is identical to spin 1/2 chains cou- 
pled by nearest-neighbor interaction J2 (J3) and next- 



nearest-neighbor interaction J\ which induces frustration 
between spins. In this model it is knownta' that when 
the ratio of the nearest-neighbor and the next-nearest- 
neighbor interaction, J2/J1 for J3 = 0, for example, is 
increased, the low energy exitation changes its behav- 
ior from the one with a finite gap like the Majumdar- 
Ghosh model (J2/ J\ = 2) to a gapless one like a Heisen- 
berg chain at Jij J\ = l/a c ;a c = 0.2411EJ (see the 
axis in Fig. 5). Similar behavior exists on increasing 
J3/J1 at J2 = 0. When J 2 = J3 (the dot-dashed line 
in Fig. 5), the model reduces to a 2D AF Heisenberg 
model with nearest-neighbor coupling J2 = J3 and next- 
nearestjifiighbor coupling J\ along only one diagonal di- 
rection.liiP Its ground state is unknown except for the 
case of Ji = J2 ^-Js , which is the simple triangular lat- 
tice with no SG£j) This is approximately similar to the 
case of the dimer model for k-(ET) 2 XE1'I3^ as shown in Fig. 
5 and the ground state experimentally observed there is 
actually an AF insulator.tf (2D Heisenberg model with 
nearest-neighbor and next-nearest-neighbctt-along both 
diagonal directions has been well studied.li3>) As seen 
in Fig. 5, 6»-(ET) 2 RbZn(SCN) 4 is located in the unex- 
plored region, yet in this region frustration is large such 
that the ground state may have a SG. This is the possible 
explanation we propose in this paper. 

Here we note that if 6»-(ET) 2 RbZn(SCN) 4 is in the re- 
gion of AF, the SP transition will be a candidate for 
the origin of the SG behavior because J\ is much larger 
than J2 and J3. It is known that the effect of interchain 
coupling on SP transition is_to enhance the tendency to- 
ward AF long range order .EZP However, the competition 
between SP transition and Neel ordering in the presence 
of frustration among the interchain couplings is also an 
interesting unanswered question. 

In SG systems, the impurity effect is interesting be- 
cause impurities easily destroy the SG and induce AF 
long range order .E3 If some kind of disorder is intro- 
duced into 0-(ET) 2 RbZn(SCN)4 without any structural 
changes, the HF calculations predict that the configura- 
tion of the magnetic ordering will be either (a) or (b) 
type AF. In recent experiments in this direction, how- 
ever, the dimeric structure was unfortunately not main- 
tained; in 6>-(ET) 2 - a; (BMDT-TTF) 2; RbZn(SCN)4 and 6- 
(ET)2Rbi_ 2; Cs 2 ;Zn(SCN)4, the structural phase transi- 
tion observed in the P-fffi. compound is supressed and 
then the SG collapses 

In summary, the origin of the SG behavior in the slowly 
cooled dimerized phase in #-(ET)2RbZn(SCN)4 is stud- 
ied both by the HF calculations and the resulting effec- 
tive dimer spin model. HF calculations indicate that the 
pairs of molecules coupled with the largest transfer inte- 
grals should be considered as dimers, and that the elec- 
trons are localized on each dimer due to the Coulomb 
interaction forming a Mott insulator. Quantum effect 
will create frustration among the localized spins which is 
the possible origin of the SG. 
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Fig. 1. Schematic donor structure of 6>-(BEDT-TTF) 2 RbZn(SCN) 4 
in the high-temperature phase (a) and the slowly cooled low- 
temperature phase (b). Each ET molecule is schematically ex- 
pressed as a rectangle. The values of transfer integrals are given 
by t = Es, where E = —10 eV and s is the value of the overlap 
integral taken from ref. 9. 



Fig. 2. Different types of AF ordering for the structure given in 
Fig. 1. Circles and ellipses represent ET molecules and dimers, 
respectively. The magnetic unit cells are shown by rectangles. 
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Fig. 3. Effective dimer model for 6»-(ET) 2 RbZn(SCN) 4 . Full cir- 
cles represent dimers and the values of the effective superex- 
change couplings are shown. Thick, thin and thin-dotted lines 
represent the spin exchange couplings, J\,J2 and J3, respec- 
tively. 



Fig. 5. The schematic phaae-i diagram of the Heisenberg spin 
model described in Fig. 3.EIP It is symmetric with respect to 
the dot-dashed line. SG, nDAFH, TLAFH, AFLRO and MG 
represent spin gap, n-dimensional antiferromagnetic Heisenberg 
model, triangular lattice antiferromagnetic Heisenberg model, 
antiferromagnetic long range order and Majumdar-Ghosh point, 
respectively. The regions where a spin gap is known to ex- 
ist are shown by the thick lines. The regions where 0- 
(ET) 2 RbZn(SCN) 4 and k-(ET) 2 X are expected to be located 
are shown by the shaded and the dotted areas, respectively. 



Fig. 4. Energy of each state, E a , shown in Fig. 2 relative to that 
of the paramagnetic solution, E p , AE(a) = E a — E p . 
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